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In this work, we examine the vahdity of quantum adiabatic theorem in thermodynamic systems. 
For a d-dimensional quantum many-body system, we show that the duration time to required 
by its ground-state adiabatic process does not depend on the microscopic details, but the scahng 
dimension of the fidehty susceptibihty da- Our resuh, therefore, provides a quantitative time scale 
of the quantum adiabatic theorem in thermodynamic systems. The quantum adiabatic theorem 
might be violated in case that the scaling dimension of the fidelity susceptibility is larger than the 
system's real dimension {da > d). 



PACS numbers: 03.65.Ca, 03.67.-a, 64.60.-i 

Quantum phase transition (QPT)jl'| is one of the most 
active research fields in condensed matter physics. For 
a quantum many-body system described by a Hamiho- 
nian H{X), a QPT occurs as its ground-state property 
undergoes a significant change at a transition point Ac. 
In order to study QPTs, people usually work on the low- 
est eigenstate of H{X). In practice, if there is no other 
mechanisms to change the lowest eigenstate, but drive 
the system from one phase to the another by changing 
the driving parameter A directly, one should ensure the 
validity of the quantum adiabatic theorem. 

The quantum adiabatic theorem states that a quan- 
tum state will not transit to the system's other states of 
different eigenenergy if the driving Hamiltonian changes 
slowly enough in time. The theorem is an extremely in- 
tuitive concept because its validity relies on the crite- 
rion of the "slowness" . This criterion, for an arbitrary 
D-level system, has been improved step by step in the 
last several decades [3, i, i, S, S 0, H ll Gol . However, 
the relation between the "slowness" and thermodynamic 
properties, such as dimensionality and various critical ex- 
ponents etc, have been paid few attention. Therefore, for 
a d-dimensional quantum thermodynamic system, how to 
define "slowness" or its relation to statistical quantities 
remains a fundamentally important question. To answer 
this question in a quantitative way is the key motivation 
of the present work. 

In this report, we start from the time-dependent 
Schrodinger equation, and show that the leading tran- 
sition probability from the ground state to excited states 
at the perturbation level is proportional to the fidelity 
susceptibility pH IT^ . which was proposed recently in 
the fidelity studies on thc QPTs 13, 14 EE [H, [13, El, 
El, im, E2, m, m, HE Si . Then we are able to use the 
scaling dimension of the fidelity susceptibility da (called 
quantum adiabatic dimension hereafter [2^ ) to quantify 
the scale of the duration time required by the quantum 
adiabatic theorem. A general inequality is established for 
the slowness criterion (in terms of the duration time) in 



the thermodynamic limit. 

We take the linear quench process, in which the driv- 
ing Hamiltonian is turned on linearly with the time t, as 
an example. In this case, the duration time tq for suffi- 
cient slowness should satisfy tq 2> kL'^" where k is inde- 
pendent of L. Therefore, if we require that a physically 
acceptable duration time is proportional to the system 
size, which is about the order of the Avogadro constant 
(6.02 X 10^^) for a realistic system, then the two limits 
of N{— L'^) oo and ro(oc N) — > oo do not commute 
with each other in case that the quantum adiabatic di- 
mension da > d (in other words oo ^ oo^ for ^ 1), 
hence the quantum adiabatic theorem might break down. 
We finally examine the validity of the quantum adia- 
batic theorem in a few of many-body systems, includ- 
ing the one-dimensional transverse-field Ising model 0], 
the Lipkin-Meshkov-Glick(LMG) model [13], and the Ki- 
taev honeycomb model |28|. In these models, da > d 
at their corresponding critical point, hence the quantum 
adiabatic theorem is violated around the critical point. 
In the Kitaev honeycomb model, moreover, the quantum 
adiabatic dimension in the gapless phase is 2-fln (for a 
size dependence of N'^ In N), which is still larger than the 
real dimension 2, so the adiabatic theorem might break 
down in the whole phase. 

To begin with, we consider a general d-dimensional 
quantum many-body system of length L and size N — 
L"^. Its Hamiltonian reads 



(1) 
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where Hj is the driving Hamiltonian, A = Ai -I- t/rQ de- 
notes its time-dependent strength with tq being the du- 
ration time scale and Xi the starting point. To be consis- 
tent with the time-dependent perturbation theory, we let 
{\4>n{t))} define the complete set of eigenstates of the in- 
stant Hamiltonian H{t), i.e. H{t)\(f>n{t)) = £n{t)\4'n{t)) ■ 
According to the quantum adiabatic theorem, the ground 
state of the system is always |0o(i)) if the driving Hamil- 
tonian Hj is turned on slowly enough (here we exclude 
those cases of the ground-state level-crossing) . Then we 
can always use the adiabatic ground state \4>oit)) to study 
QPTs in the parameter space of A. 
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According to quantum mechanics, the system's state 
can be expressed as a hnear combination of the adiabatic 
eigenstates, 



(2) 



which is required to satisfy the time-dependent 
Schrodinger equation, i.e. 



^^^\^{t))=H{t)\^it)). 



(3) 



Here we set h = I. Introducing the unitary transforma- 
tion 



(4) 



an{t) = an{t) exp -i / e„t dt 



and combining Eq. ^ and Eq. ^ together, we can 
obtain 



However, it is not easy to estimate exactly the vahies 
of the integrals in Eqs. (O and ([7]). To see the quahta- 
tive behavior of the leading term of the fidehty, we first 
make use of |0„(A)) as reference states, then the energy 
levels vary slowly with time. Under this approximation, 
the fidelity is the same as the perturbative form of the 
Loschmidt echo 01 
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Here the second term denotes the transition probability 
and a phase factor from Eq. Q has been normalized out. 

The second alternative approach is to find the bound 
of the integral in Eq. ([7]). Because of 



exp 



1, 



(10) 



we have 



dcLr. 



dt 



d„i((t)m\dt<j)m} 

E{4>rn\dtH\(j}n)an 
exp 
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where dx4> = d(j)/dx and ujnm = Cn — Cm- The quan- 
tum adiabatic theorem is based on the approximation 
that if the second term on the right hand side of the 
above equation is small enough compared with the first 
term, then the mth state will keep its position except 
for an accumulation of a Berry phase from the first term 
of Berry connection. Such an approximation holds true 
for a finite-size system with a finite-size energy gap, but 
should be treated very carefully for a thermodynamic 
system in which the gap might vanish and long-range 
correlations appear in various distinct ways. 

Now we suppose the Hamiltonian evolves from A to 
X + SX during a finite time interval At, that is At = 
SXtq. The SX is small enough for the validity of the time- 
dependent perturbation theory. We will return to this 
requirement later. At time t — 0, — I, a„i — 0, so the 
system is at the ground state (j)o{t ~ 0), then at t = At, 
we have, to the first order, 
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where i/p™ 



{(t)n\dxH\(l)m)- To see the vahdity of 
the adiabatic theorem, we need to address the fidelity 
between \4>o{t)) and |5'(t)}- For the normalized states 
\(j}o{t)) and |*(t)}, the fidelity is 



The adiabatic theorem requires that F : 



(8) 
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Then we obtain a lower bound of the fidelity 
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hence an upper bound of the transition probability. 

Mathematically, Fi defines a distance between |0o(^)) 
and Y^{t)) and F2 a distance between |0o(A)) and |(/)o(A-|- 
6X)). Therefore, Fi, F2, and F in Eq. © for a "triangle" 
in the parameter space. Our concern here is that the 
transition probabilities in both f 1 and F2 are determined 
by the fidelity susceptibility [HI, 
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which defines also the scale of the original fidelity F de- 
fined in Eq. ([5]). In previous studies on the quantum 
adiabatic theorem, the formulism given in Eqs. (|S] |12p 
are familiar to us, however, few attention has been paid 
to the scaling behavior of the quantity (the fidelity sus- 
ceptibility) until recently [H, [l^ . 

For a c?-dimensional system, the fidelity susceptibility 
of the driving Hamiltonian has its own dimension da [23| 
instead of the system's real dimension though in many 
cases both dimensions are equal. That is 



XF cc L 



(14) 



given that L is larger enough. In the critical region, the 
quantum adiabatic dimension da = 2d + 2C — 2Av[l6j 
with d, C, and Ay being the real dimension, dynamic ex- 
ponent and scaling dimension of the driving Hamiltonian 
respectively. Clearly, in this case, the quantum adiabatic 
dimension da can be larger than d. In the non-critical 
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TABLE I: The adiabatic dimensions d^,(F^, and ( above, 
at, and below the critical point) for the ID Ising model 
[14 ]. the LMG model [S^l, and the Kitaev honeycomb model 
(KHM) [2l|,[2|. 



Model (critical point) 


d 








ID Ising model [he = 1) 


1 


2 


1 


1 


LMG model(/!,e = 1) 


1 


4/3 





1 


KHM(Jc = 1/2) 


2 


5/2 


2 


2+ln 



region, the correlation length is finite, then we usually 
have da = d or da < d. For instance, in the fully po- 
larized phase of the LMG model, da = [13 (here the 
LMG model is considered as a one-dimensional system 
with infinite-range interactions). In Table HI we show the 
adiabatic dimension for three exactly solvable models, i.e. 
the one-dimensional transverse-field Ising model the 
LMG model |27| , and the Kitaev honeycomb model [1^ , 
around their corresponding phase transition point. These 
data are collected from the recent fidelity approaches to 
QPTs, as shown in the caption of the table. 

However, the results obtained from the perturbation 
theory is valid only if the change in the driving Hamil- 
tonian is very small. Physicists are interested in the 
case that the change in the Hamiltonian, though varies 
slowly with the time, is not small, then the perturbation 
method can not be applied directly. To solve this prob- 
lem, we first assume, without loss of generality, that the 
system evolves from Aj to Xf during the time interval 
To, and the quantum adiabatic dimension in this region 
is da- Secondly, we divide the interval A/ — A.; into M 
subintervals (Here M defines also the scale of tq), that is 
6X = {Xf — Xi)/M . The total leading transition proba- 
bility to excited states Pt then scales like 



XF- 



Because of xf oc L"*" , M should be at least about the 
order of L''" to ensure the validity of the perturbation 
method. The final (minimum) probability of staying in 
the ground state \(j)o) at A/ then becomes 
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(15) 



The quantum adiabatic theorem requires ~ 1. In the 
thermodynamic limit, we arrive at a simple inequality. 



To » kL" 



(16) 



where k is a L-independent quantity. The above inequal- 
ity concludes the key result of the present work. 

Clearly, the inequality (fTBl) defines the scale of dura- 
tion time required by the quantum adiabatic theorem. 



That is, a sufficient adiabatic condition should defines 
the duration time in order of L"*" . For a realistic (ther- 
modynamic) system, the only comparable scale is the 
system size N = L"^, which is about the order of the 
Avogadro constant (6.02 x 10^"^). If we require that a 
physically acceptable duration time is, at most, propor- 
tional to the system size, then the two limits of N — > 00 
and To — > 00 do not commute with each other if da > d, 
and the quantum adiabatic theorem might be violated. 
On the other hand, in case that the starting and ending 
points are located in two regions with different quantum 
adiabatic dimensions, say dai and da2, the required du- 
ration time is determined by tq ^ ^^max(d„i,d„2)^ 

Now we take the one-dimensional transverse-field Ising 
model as an example to illustrate the validity of the quan- 
tum adiabatic theorem according to our criterion. The 
Hamiltonian of the Ising model reads 
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-E 



(17) 
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where h{t < 0) = —t/rQ. The Hamiltonian p7|) has 
been used as a prototype model in both studies on the 
ground-state fidelity [l^ and dynamics of QPTs [s^, [3l| . 
If the quantum adiabatic theorem holds true, t < —tq 
corresponds to the paramagnetic phase, and — tq < t < 
is the ferromagnetic phase. A second order QPT occurs 
at t = —Tq. The fidelity susceptibility in the whole region 
can be calculated as 



XF 



E 

fe>0 



d9k 
dh 



with k^ir/N, Sir/N, ■ ■ 

Mk 1 
dh 



,n{N -1)/N, and 
sin k 



2 1 + /i2 - 2ft, cos fc ■ 



(19) 



(20) 



It can be shown that xf oc N"^ for h = 1, while xf oc 
N for h ^ 1. Therefore, according to our criterion, if 
the starting and ending point are in the same phase, i.e. 
^ii^f) < ""^0 or —To < ti{tf) < 0, the duration time 
required by the adiabatic condition is N <^ tq. However, 
if ti < — Tq and — tq < tf < 0, the system will across the 
transition point h = 1, at which the adiabatic dimension 
is 2. So the required duration time should satisfy tq » 
kN^. This observation is consistent with result obtained 
via the Landau-Zener formula 0, 13| in the recent studies 
on quench dynamics in the Ising model [sij . Therefore, 
the quantum adiabatic theorem might break down at the 
critical point. 

However, the problem is still subtle because the transi- 
tion point is not a region but a "point" . Then the large N 
behavior might be quite different from that of the infinite 
limit. In Fig. [U we show the scaling behavior of the fi- 
delity susceptibility around the critical point. We can see 
that only at the critical point, xf/N cx N. While away 
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FIG. 1: The scaling behavior of the fidelity susceptibility 
around the critical point of the one-dimensional transverse- 
field Ising model. 

from the critical point, though the closer to the critical 
point, the larger the fidelity susceptibility, the later will 
be finally saturated to 

X£ ^ f for /i< 1 

respectively as N increases. Therefore, for any simulation 
on a large but finite sample, the duration time should 



satisfy tq ^ kN"^ in the region close enough to the critical 
point. While in the infinite N limit, the condition tq ^ 
kN'^ is valid rigorously only at the critical point. 

Moreover, we can see from Table I that the quantum 
adiabatic theorem does not hold true at the critical point 
of both the LMG model and Kitaev honeycomb model 
also. On the other hand, it has been found recently that 
the quantum adiabatic dimension in the gapless phase of 
the Kitaev honeycomb model is 2+ln 2^ , which is larger 
than the real dimension 2. Our quantum adiabatic con- 
dition implies that the quantum adiabatic theorem might 
be violated in the whole gapless phase of the Kitaev hon- 
eycomb model. 

In summary, we have proposed the quantum adiabatic 
condition for quantum systems in the thermodynamic 
limit. A general inequality between duration time re- 
quired by the quantum adiabatic theorem, the system 
size, and quantum adiabatic dimension are established. 
For the commonly studied linear quenches, our results 
show that the adiabatic condition might be violated if 
the adiabatic dimension is larger than the real dimen- 
sion. This phenomenon usually occurs at the quantum 
critical point and those strange phases of da > d. 
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